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ABSTRACT 

The spectral problem of thin elastic shells in membrane approximation does not satisfy the 
classical properties of compactness and so there exists an essential spectrum. In the first 
part, we propose to determinate this spectrum and the weakness directions in the shell. 
We particularly study the case of homogeneous and isotropic shells with some examples. 
In the second part, we consider an elementary model problem to study the propagation of 
singularities and their reflections at the boundary of the domain. In the last, we study the 
problem of propagation for an isotropic cylindrical shell and we show that the equation of 
propagation does not depend on the Poisson coefficient. 
keywords : shell, essential spectrum, propagation of singularities 



1 Introduction 

1.1 Classical and non classical vibrating problems in shell theory 

We consider a thin shell with a middle surface S. This surface is described by the map, 

y=(y 1 ,y 2 )en^r(y 1 ,y 2 )eR i (1) 

where Q is a domain of the plane. Let u(y) = (ui(y) , u 2 (y) , u 3 (y)) be the displacement vector 
of the surface and its covariant components. We introduce the Hilbert space H = (L(Q)) 3 
and we denote by (u,v) the scalar product. The displacement u belongs to the subset 
V\ = H 1 ^) x H l (Q) x H 2 (Q) which can be modified to take boundary conditions in account. 
The variational form of the problem of vibrations is (cf.[12, ch. VI]), 
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We search for u e V\ so that, 

Vi> e Vi, a m (u,v) + e 2 af(u,v) = X(u,v) (2) 

The bilinear forms a m and ay correspond respectively to the membrane problem and the 
flexion problem. They are continuous on Vi. We denote by A the spectral parameter. This 
problem is classical with a selfadjoint operator and compact resolvent and so there exists a 
sequence of eigenvalues 

O < A < Xi < ... < X k < — ► +00 (3) 
with orthogonal modes (cf.[12, ch.I]). 

If the relative thickness of the shell, e, is very small, then the membrane approximation is 
an appropriate representation. The formulation of this problem is different. In this case, u 
belongs to the space V = x x L 2 {Vt). The inclusion of V in H is dense and 

continuous but is not compact. The problem is written as, 

We search for u £ V so that, 

Vf G V, a m (u,v) = X(u,v) (4) 

This spectral problem is an elliptic system with mixed order. The classical properties of 
compactness are not satisfied and the spectrum both contains a sequence of eigenvalues 
depending on the domain, and an essential spectrum. 



1.2 Essential spectrum 

Let if be a Hilbert space and A a selfadjoint operator. The resolvent set is defined by, 

p(A) = {(/ (A-CIdy'eCiH)} (5) 
Its complement, the spectrum is constituted, 

- of isolated eigenvalues of finite multiplicity; for these £, (A — (Id)" 1 does not exist. 

- of other values for which (A — (Id)" 1 exists but does not belong to C(H). They 
are eigenvalues of infinite multiplicity, accumulation points of eigenvalues and continuous 
spectrum. 

The set of these ( which are not isolated eigenvalues of finite multiplicity is the essential 
spectrum S ess (A). 

It can be characterized as the set of ( for which there exists a sequence (uk) called Weyl's 
sequence so that, 

1 1 u k 1 1 =1 

Uk — > in H weakly (6) 
(A — (Id)(u k ) — >0 in H strongly 



For very small data, we can obtain a large response and so these sequences can be physi- 
cally interpreted as some kind of resonance. This local phenomena are quick oscillations in 
some directions which are called weakness directions. We will see that the singularities will 
propagate along directions which are orthogonal to these weakness ones. 

In the first part of this paper, we study the essential spectrum of a vibrating shell in mem- 
brane approximation. We show how it is possible to determine this set and the correspond- 
ing weakness directions and we give some examples. In the second part, we introduce a non 
classical model problem and its essential spectrum and we investigate the propagation and 
reflection of singularities. Finally, in the last part, we study these problems for a vibrating 
cylindrical shell. 



2 The case of shells in membrane approximation 

Let S be the middle surface of a thin shell, described by the map, 

(yVle^^VjeR 3 (7) 

We define the fundamental forms, 

A(x\x 2 ) =a a/3 x a x 13 (8) 

and 

B(x\x 2 ) =b af3 x a x (3 (9) 

The equations of the vibrating shell in the membrane approximation give the following 
spectral problem, 

-D l T n -D 2 T 21 =Xu 1 

—DiT 12 — D 2 T 22 =\u 2 (10) 
-b aP T<*P = \u 3 

where u 1 , u 2 and w 3 are the contravariant components of the displacement and T a/3 , those 
of the stress tensor (cf.[ll, ch. X], [1]). The covariant derivatives of tensor T af3 are given by, 

using the Christoffel symbols. We define the strains of the shell by coefficients 7a m («), 

7a m H = t^(D x u^ + D^ux) - b x ^u 3 (12) 

where D\ denotes covariant derivative on (S) and the behaviour of the shell is obtained with 
the elasticity coefficients a Q/3AM so that, 

T afS = a al3X ^ x ,(u) (13) 
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We will note A = [a a/3AAl ] the stiffness matrix. 

By replacing the expressions (|13J1 of T a/3 in (|1(J|). we obtain an explicit spectral problem 
on the displacement u. There appears derivatives of second order in u\ and w 2 and of first 
order in u 3 . The classical properties of compactness are not satisfied and there exists an 
essential spectrum. We have weakness directions noted by (£i,£ 2 ) and the orthogonal ones 
(xi, x 2 ) will be the directions of propagation of the singularities. The values of A and (£1, £ 2 ) 
correspond to the non-ellipticity of the system in the Doughs and Nirenberg sense (cf. [11 
ch. Ill] [5]). In this spectral problem, we have three equations corresponding to (JIUj) . The 
equation i is of order (sj + Sj) in the variable Uj with si = s 2 = 1 and s 3 = 0. 
To determine A and (£i,£ 2 ), we write that they are solutions of the equation we obtain by 
writing that the determinant of the principal symbol is zero: 

/ a^Ud a al ^U(3 ia al ^b^« \ 
det a^Up a a2 P 2 Up ta a2 ^b^ a = (14) 
\za al ^b c ^ a ia^b^a X-a a ^b aP b c J 

which we write in the condensated form, 

A 12 A 13 \ 

= (15) 





(A xx 


A 12 


A 13 


det 


A 12 


A22 


A 23 




Ul3 


A 23 


A- B; 
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Whe then have 



> „ , A 2 13 A 22 + A 2 23 A U - 2A 12 A 23 A 13 

33 + A~A~ 42 U°J 

^11^22 — ^12 

The numerator B 33 (A U A 2 2 — Af 2 ) + A 2 3 A 2 2 + Al 3 A n — 2A 12 A 23 Ai 3 is polynomial of degree 4, 
homogeneous in £i,£ 2 . By calculating every coefficient, we see that all of them are multiple 
of detA and finally that it is exacty equal to detA . [B(— £ 2 , £i)] 2 . It is easy to calculate the 
denominator, 

A n A 22 - A\ 2 = c 22 tf + c n £ 2 4 + (C33 + 2c 12 )£ 2 £ 2 - 2C23&2 - 2c 13 Ul (17) 

where c a /3 are the cofactors of matrix A. We then obtain the relation between A and (£1, £ 2 ). 
By replacing the components of a vector of the weakness direction by those of the direction 
of propagation (xi,x 2 ), we have, 

{X!,X 2 ) = (18) 

and then, 

A = '^""'l' (19) 

S^^XaXpXxXfj, 

where s a/3Xfl are the coefficients of A -1 (compliance coefficients). 

Let us recall that the coefficients which appear in that last expression depend on (y 1 ,?/ 2 ). If 
a point is given on S, then the spectral parameter A belongs to a segment 

S ess =[A l5 A 2 ] (20) 
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and the whole essential spectrum is the set of all these segments when (y 1 ,!/ 2 ) draw fl. 
By noting that A is a quotient of two quadratic forms in X = (x\ , x 2 , Xix 2 ), 

A = (21) 

l XQ 2 X 1 ; 

we can write that, 

£ ess = [Ai, A 2 ] C [min (eigenvalue of Q2 1 Qi) ) max(eigenvalue of Q^Qi)] (22) 

and we obtain the classical inclusion (cf.fll ch. XI]), 

Sess C [0, B 33 ] = [0, a aK % ap b Cri ] (23) 

We note that is reached in every hyperbolic point of the shell but not in elliptic ones as it 
is easily seen from (jl9|) . Moreover, there are some cases for which B 33 is not reached. For 
instance, if we consider an isotropic cylindrical shell (x 1 = Ry 1 ; x 2 = Rcosy 2 ; x 3 = Rsiny 2 ), 
then we have 

E - = [0, * [0, B 33 = R2{1 E _ u2) ] (24) 

Let us now consider the case of an isotropic shell. The elasticity coefficients are given by, 

E , . o„ „„, fli 2v 



and so we have, 



and 



a a ^ = 2 (a aA a^ + a^a^ + a Q V^) (25) 



detA = f (a*V 2 - (a 12 ) 2 ) (26) 

2(1 - u 2 )(l + v) 

s^^XaXpXxXp = \- {ax\x\ + a 22 x\ + 2a\ 2 x\_x 2 ) 2 (27) 
hi 



where [a a p\ = [a af3 ] 1 . Finally, we obtain the following outstanding form of (|19|). 

B(X 1 ,X 2 ) T 2 



A = E 



A(x 1 ,x 2 ] 



(2i 



where A et B are the two fundamental forms of the surface S. It appears that the essential 
spectrum depends only on the geometry and the Young modulus but is independent of the 
Poisson coefficient. A geometrical interpretation of the quotient of the fundamental forms 
MgliggJ j s the normal curvature k x of the surface in direction (x±, x 2 ). The essential spectrum 
is then exactly the segment, 

S ess = [ E.Infkl E.Supkl ] (29) 

Conversely, if A G T, ess is given, then we can find the couples (xi,x 2 ) which define the 
directions of propagation of singularities. 

In an elliptic point of S we have two directions (xi,x 2 ) but in a hyperbolic point several 
cases are possible 
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Let us note, k\ = Infk x < < k 2 = Supk x which correspond to the principal curvatives. 
We have the following results: 

1. If we suppose that —k\ < k 2 then, 

If A e] 0, Ek\ \, then there are four directions of propagation 
If \ e] Ek\,Ek\ [, then we have only two directions, 
and if A > Ek\, then there are no direction. 

2. If —k\ > k 2 , then there are two directions of propagation if X e] 0,Ek 2 [ and zero if 
\>Ek\. 

Example: 

We consider the hyperbolical paraboloid shell defined by the map, 

(y\ y 2 ) — (x 1 = y\ x 2 = y\ * 3 = | Uf - (y 1 ) 2 ]) (30) 



we easily calculate, 



an = 1 + ^ (y 1 ) 2 ; ai 2 = -^J/V 5 «22 = 1 + ^ (y 2 ) 2 (31) 



&n = -&22 = - - , ; &12 = (32) 

o 2 V°n a 22 - af 2 



Denoting by ^ the polar angle of the weakness direction £, ip e] — f , f [, we obtain 

L A(tanip, 1) 
For 6 = 0, 5, c = 0, 1 and y 1 = y 2 = 0, we have 

-0,4(tanV> - 1) 



(33) 



X = E 



tan 2 ip + 1 



(34) 



and then S ess = [0, 0.16]. In this case k\ = —k 2 and we have four directions (which could 
coincide) which are symmetrical about the polar axis as it is shown on the following figures. 
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For A = 0, we have two double directions. This case corresponds to the static problem and 
the directions of propagation (which are the same that the weakness ones), are also those of 
the asymptotic curves of the surface (bu = — 6 2 2 ; &12 = 0), 




Then when A increases, we obtain four directions, for example for A = 0, 04. E, 




and for A = 0, 12.E, 




When A reaches the upper bound of the essential spectrum, A = 0, 16. E, two directions 
disappear and the others become coincident: 



Let us consider another point y 1 = 2, 5 and y 2 = 5. We have, 



A = E 



0,2(tan 2 ip - 1) 
3(5tan 2 ^ - Atamp + 2) 



(35) 



and S ess = E [0, 2.36 1CT 3 ]. Here, fcj = 0.232 lO" 3 and if A G [0, 0.232 10- 3 £[, we have four 
directions: 

For A = then ip = 45° and we have the two double directions of the asymptotic curves as 
in the previous case, 




then we obtain four directions. For A = 0.1 10~ 3 E then = -75,2°, ip 2 = -26,8°, 
^ 3 = 39, 6°, ^ 4 = 54°, 
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For A = 0.232 10 3 E = Ekf, two directions become coincident in ip 1 = t/> 4 = 73° and 
^ 2 = -18,4°, ^3 = 37,2°, 




For the values of A larger than 0.232 10 3 E, there are only two directions: 
If A = 0.232 10" 3 £ + then = ^4 disappear, 




for A = 2 1Q- 3 E then ip 2 = 10, 2°, ip 3 = 23, 6° 




and finally for A = 2.36 10 3 E the two directions become coincident in ip 2 — i>3 — ^7,4°, 
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and then they disappear because we go out of the essential spectrum. 



3 Propagation of singularities in a model problem 

The equations of vibrating elastic thin shells in membrane approximation are rather com- 
plicated. So to give a good idea about the properties of propagation of singularities we will 
first study a model problem (cf.[13]) and then we will look at the case of a particular shell. 
We consider the following spectral problem, 

f -Ami + bu 2 ,2 -Xui = fi ^ 
\ -bu 1:2 + cu 2 -\u 2 =h 



where the unknowns (ui,u 2 ) are functions of the two variables y = {y l -,y 2 ) in a domain Vt 
of the plane. We suppose that b and c are given real numbers and f — (fi, f 2 ) is a right- 
hand side which will be defined later. The boundary conditions are for example Dirichlet 
conditions U\ = on the boundary dQ. 



3.1 Essential spectrum 

We define the spaces H = L 2 {Vl) x L 2 {Vl) and V = Hq(Q) x L 2 {Vt) and we can write the 
problem in the form, 

(A -\)u = f (37) 

where A is a selfadjoint operator in H. 

The problem involves derivatives of second order in ui and of first in u 2 . This problem of 
mixed order does not satisfy the classical properties of compactness (obviously the inclusion 
of V in H is not compact). We determine the essential spectrum by writing that this 
problem is not elliptic in the sense of Doughs and Nirenberg. We have two equations and 
two unknowns and the indices are s± — 1 and s 2 = (cf.fll ch. Ill], [5]). We look for the 
existence of a nonzero real vector £ 2 ) so that the determinant of the principal symbol 
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det ( $ \f 2 %K \ ) = (38) 



vanishes, 

A P + ( £1 

-ife£ 2 c — A 

and it is easily checked that A must belong to [c — b 2 ,c\. That segment is the essential 
spectrum which does not depend on the domain Q. 

In general, the set of eigenvalues of a selfadjoint operator in a Hilbert separable space is denu- 
merable and the corresponding eigenvectors are orthogonal. So, in particular, the eigenvalues 
of A contained in the essential spectrum have measure zero. 

Let us consider for instance that Q is the square [— | , | ] 2 , with Dirichlet conditions on the 
edge. The functions, 

Ul (y\y 2 ) = A lC os{2p + l)y 1 .cos{2n + l)y 2 (39) 

u 2 {y\ y 2 ) = A 2 cos{2p + l)y\sin{2n + l)y 2 (40) 
are eigenfunctions of the problem when the following equation is satisfied: 

« 2 - 1 + i^±ii!_^ (41) 



c-A (2n + l) 2 (2n + l) : 

with 

b{2n + l)A x + (c - X)A 2 = 

Then, when the integers p and n are given, we find two eigenvalues which are at the inter- 
sections of a straight line and a hyperbola (see (j41|) V One of them is in ]c — b 2 , c[ and the 
other is larger than c. For any p and n we can show that the subset of these eigenvalues is 
dense in the essential spectrum and also in ]c, +oo[. 



3.2 Propagation of singularities 

For A belonging to }c — b 2 ,c[, there exist two different directions (£1,^2) which satisfy to (j3HJ). 
The orthogonal directions are the two directions of propagation (cos9, sinO) with, 

tan 2 6 = A ~ (C ~ &2) (42) 
c — A 

At the extremities of the essential spectrum (A = c — b 2 , resp. A = c), the two directions of 
propagation become coincident {6 = 0, resp. 6 = |). 

We can note that the substitution of 1*2,2 into the first equation gives the wave equation, 

A — (c — b 2 ) b 

- ux n H Mi,22 - Ami = fi -/ 2 ,2 (43) 

c — A c — A 

whose characteristic directions correspond to the directions of propagation which have been 
defined by 9 (cf. [4]). 
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If A is in the essential spectrum but is not an eigenvalue of A, then the range of A — A is 
dense in V = H^ 1 ^) x L 2 (Q) and A — A is injective. For some right-hand sides /, there 
exist solutions which are unique but the resolvent is not continuous in these spaces. 
Let us consider problem with the right-hand side f\ — 0; f 2 — 6(y 1 )6(y 2 ). This right- 
hand side does not belong to V but we can consider this problem as the research of a 
fundamental solution (cf.[6]). 

As A is not an eigenvalue, if there exists a solution, then it is unique. 
When A is exterior to the essential spectrum, equation (J43)) becomes 

- «i,u + A ~ (C ~ 6 \ ,22 - A« a = -^-5{y l )5\y 2 ) (44) 
c — A c — A 

and is elliptic and the right-hand side belongs to H S (Q) with (s < —2) (cf.[10]). If the 
solution u\ exists, then it belongs to H S+2 (Q) (and u 2 to H S+1 (Q)). In that case there is no 
propagation. 

Let us suppose now that A belongs to ]c — b 2 , c[ and is not an eigenvalue, problem (J35J) can 
be also written in the form, 



-Ami + bu 2 ,2 -Ami = , , 
-bu lt2 + cu 2 -\u 2 = S(y 1 )5(y 2 - yHanQ) 

where tand is the slope of one of the characteristics {Dq). We will look for solutions in 

the form of an asymptotic expansion of the singularities across the characteristics. To take 
orders of derivatives in consideration, we propose the following expansions, 

Ul (y\ y 2 ) = Uliy^iy 2 - yHanQ) + Ul{y x )Y{f - yHanQ) + ... . . 

u 2 (y\ y 2 ) = UlW)8'{y 2 - yHanQ) + U^y^y 2 - yHanQ) + ... l40j 



where Y is the Heavyside function. We substitute these expressions in the problem and we 
can identify the leading terms, in 5"(y 2 — y 1 tan9) in the first equation and S'(y 2 — y 1 tan8) 
in the second. 
We obtain a linear system, 

'l + tan 2 6)U\ -bUl =0 , . 

-bU\ + {c-\)Ul =0 { ' 



which admits nonzero solutions because the determinant vanishes when A belongs to the 

essential spectrum and Q is defined by fl4^j) . 

The identification of the next terms gives another linear system, 

[l + tan'Qpfiy^-bUKy 1 ) =2tan6 d ^(y 1 ) 

_ bU 2 {y l ) + {c _ X)U 2 {y l ) =S{y l ) 
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with the same determinant. To obtain solutions, the right-hand sides have to satisfy the 
following compatibility condition, 



dU 1 

2{c-\)tan9 —\{y 1 )+b5{y 1 ) = Q (49) 
ay 

and then we can find function U\ which defines the propagation of the singularity along the 
characteristic (Dg) 

U ^ = 2(^k^ {Ylyl)+C) <50) 
where C is arbitrary and can be obtained by the Dirichlet conditions on the edge of the 
domain. 

In that case, we can find a solution u\ of (JHfij) belonging to if s+ 2(f2), [s < —2). Then the 
solution Ui belongs to H S+ 2(Q). The singularities of these solutions are then much more 
important than in the case when A is out of the essential spectrum and moreover we saw 
that this singularity can propagate. 



3.3 Reflection of the singularities 

The propagation of singularities is characterized by U\ and is along the characteristic straight 
line which cuts the edges of the domain in two points. So we will have two conditions to 
determine only one constant C. It shows that the singularity will not disappear by reaching 
a point P on dQ. We will have a reflection on another characteristic the slope of which is 
—tanO. 




t 



Let (P, t, ft) be a local referential. The vector t is tangent to dfl and n is normal and exterior. 
We denote by 4> the polar angle of fx in the plane (y 1 , y 2 ) and (n,t) are the coordinates of a 
point in this local referential. By shifting the origin of y in P, we have y l = tsincf) + ncoscf) 
and y 2 = —tcostfi + nsincj). We write that the displacement is the superposition of two 
singularities, incident one U\ and reflected one u±, which gives the leading terms on the form 

Ul{y l )5{y 2 - yHanO) + U^My 2 + yHanO) (51) 
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By using the new coordinates, we obtain on the boundary dtt, 

UKtsincj)) U{(tsin(f)) 



(- 



+ 



-cos(j) — tanO.sincj) —coscj) + tan9.sin(f> 
which must vanish. Therefore, 

U}(tsin(j)) U\{tsin(f>) 



)S(t) 



cos((f) -9) + cos((f) + 9) 







(52) 



(53) 



which shows that if the incident characteristic is almost tangential to the contour of the 
domain, then cos(0 — 9) is small and the intensity of the reflected singularity is very strong. 
In certain cases, the mechanism of propagation and reflection may improve the intensity of 
the singularities. 

By introducing the expressions of U\ and U\, we easily obtain the following relation, 

Ccos((f) + 9)- dcos((f) -9) = 2sin<p.sin9 (54) 

Obviously the reflected singularity on the characteristic {D_q) will be reflected again when 
this straight line will reach the edge in another point and so on. 

As a first example, let us suppose that Vl is a disk and that 9 — |. By following the 
characteristics from a point inside Q, we see that a singularity will be propagated and 
reflected along the sides of a rectangle which is inscribed in the circle dil. 




We denote by Ci, C 2 , C 3 and C 4 , the constants corresponding to the four sides and we have 
the following system, 



C lC os(0 + f ) - C 2 cos(0 - f ) 
C 2 co S (0-f) + C 3 co S (0 + f) 
C 3 co S (0+f)-C 4 co S (0-f) 
C 4 co S (0-f)-C lC o S (0+f) 



V2sin(J) 
-\/2cos4 
\f2sin(f) 
-\plcosi, 



(55) 



The determinant is zero and the compatibility condition is satisfied. 
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This shows that the asymptotic structure (or at least its leading terms) exists and is not 
unique. This suggests that the corresponding value A = c — y is an eigenvalue. It is 
certainly the case for b 2 = 2c because A is equal to zero and the homogeneous equation 
which is associated to admits some nonzero solutions vanishing on dQ. 

Let us consider now a second example in which Q is the square [— | , | ] 2 . We may have a lot 
of successive reflections on the edges and we obtain an angled line. 



We denote by C\, C 2 , ... Ck, ... the constants and we have the following system for all value 
of 9 different from multiples of ~ , 



C\ — C*2 


= 


C*2 + C*3 


= -2 


C3 — C4 


= 


C4 + C5 


= -2 


etc... 





V 



It is not hard to see that the piecewise straight characteristic is closed if and only if tanO is 
a rational number. In this case we obtain a linear system the determinant of which is zero 
as in the previous example and the compatibility condition is satisfied. 

For other value of 9 (or other Q) the corresponding trajectory of the reflected singularities 
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may be somewhat entangled, so implying complicated phenomena of resonance. 

Let us remark finally that at the extremities of the essential spectrum, (A = c—b 2 and A = c), 
the two directions of propagation are parallel to the axes of coordinates. The incident and 
reflected directions would be the same so the reflection does not make sense. 



4 The case of the cylindrical thin shell 

We consider the cylindrical shell defined by the map, 

(y 1 , y 2 ) E [0, 1] x [0, 2tt[ — ► (x 1 = Ry\ x 2 = Rcosy 2 , x 3 = Rsiny 2 ) (57) 

We calculate, 

an = a 22 = R 2 (58) 

b 22 = R (59) 

The other coefficients of the fundamental forms and the Christoffel coefficients are equal to 
zero. 

The essential spectrum is the segment [0, Jj-]. Denoting by A a value in the essential spectrum 
and by 6 the polar angle of the directions of propagation x, we obtain 



A _ E 



tnn 2 i 2 

(60) 



.1 + tan 2 6 

The equations of the vibrating shell in the membrane approximation are, 



T^+T 12 ^ =-Au 1 

T 12 , 1+ T 22 , 2 =-\u 2 (61) 
RT 22 =-Am 3 

where u 1 , u 2 and u 3 are the contravariant components of the displacement and T alS , those 
of the stress tensor. By using the stiffness matrix and the strains 7a m we have, 

T 11 = K(u hl + vu 2 , 2 - uRu 3 ) 

T 22 =K(uu 1A +u 2>2 -Ru 3 ) (62) 

T 12 =tf^(ui,2 + U2,l) 



where K = j^^. 
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We obtain the following homogeneous spectral problem, 



Mi,H + ^""1,22 + ^^2,12 - vRu 3j i 

^Ml,12 + ^^2,11 + U 2 ,22 - Ru 3 ,2 



(63) 



To study the propagation of the singularities as in section 3, we introduce the following 
right-hand side which represents a point normal force 

"1,11 + ^"1,22 + ^«2,i2 - vRu 3>1 + j&pUi = 

^1,12 + V M 2,H + "2,22 - RU 3 ,2 + THpUi =0 (64) 

uu lfl + it 2 ,2 - Ru 3 + wrU 3 = 1 ^5(y 1 )5(y 2 ) 



We suppose that A belongs to the essential spectrum but is not an eigenvalue. We search for 
asymptotic expansions of the displacements on the form, 

u a (y\y 2 ) = U 1 a (y 1 )6(y 2 -my 1 ) + U 2 a (y 1 )Y(y 2 -my 1 ) + ... a = 1,2 
u 3 (y\y 2 ) = t/ 3 V W - my') + U 2 {y')8{y 2 - my 1 ) + ... ^ 



where m = tan6 is the slope of the direction of propagation, and we substitute these expres- 
sions in the problem. By identifications of the leading terms, we obtain the following linear 

system, 



2 



-m 



^UKy^-m^U^ + umRU^y 1 ) =0 
- W) + (¥™ 2 + WV) - ^sV) =0 (66) 



-vmUKy^ + UKy^ + ^-R^liy 1 ) =0 



We shall have nonzero solutions if the determinant of the matrix of this homogeneous system 
vanishes. We then obtain a relation between m and A, 

9/ 9n m 4 Ehm 4 , , 

A = KR 2 1 - v 2 ) = 67 

(l+m z ) z R z (l + m z ) z 



which is the same as ()60|) which was given by the Doughs and Nirenberg method. We note 
that it does not depend on the Poisson coefficient. It is easy to calculate U\{y l ) and U^y 1 ) 
according to U^y 1 ): 

UTf}^ 1 — 1 

= -^TJ^-^V) (68) 

and 



U 2 {y)-R (1 + m2)2 U 3 (y) 



(69) 
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Let us write the identifications of the next order terms, 

+ V)^! (V 1 ) - ^UKy 1 ) + umRUKy 1 ) = 2mU\ V) - x -^-XJ\ '(y 1 ) + uRU^ '(y 1 
±±^V) + (±f±m 2 + 1)[/ 2 V) - RUl(y l ) = +(1 - v)m\J\ V) - V) 
-umUt^) + f/ 2 V) + - R)U2(yi) = -vU{ V) + j^5(y^ 

(70) 

The determinant of the system is zero and then we have to satisfy a compatibility condition. 
Denoting by the cofactors of the matrix of the linear system and 6, the three right-hand 
sides of (J7UJ), this condition is, 

61M13 - & 2 M 23 + 63M33 = (71) 



with the same determinant. A staightforward computation gives the cofactors, 

Mi 3 =^(z/m 3 -m) 

M 23 =^((2 + z/)m 2 + l) (72) 
M 33 = ±^(m 2 + l) 2 



and then the right-hand sides have to satisfy the following compatibility condition, 
vmR{vm 2 - V)U\ \y l ) + — [{v 2 +v- 4)m 2 +3u- 1)U% '(y 1 ) + 

(vrrt - (1 + " ) 2 (2 ~" ) m 2 + l -^-)ul V) = -(m 2 + l) 2 ^*/ 1 ) (73) 



We replace U^(y l ) and U\(y l ) by their values (cf. (|68p and ()69jl) . and we finally obtain 



1 

'in- I ' " ' "' ~ 7 



777 

4i2(l - z^ 2 )^^^ 1 V) = K + l) 2 ^^?/ 1 ) (74) 



and then we can find the equation of propagation. By replacing the value of K (cf. (J52Jl ). we 
have, 

u ' = ^Te^ 1 ^ (75) 

and we see that this equation does not depend on the Poisson coefficient. 
The function U$ which defines the propagation of the singularity is 

^V) = ^^0V)+C) (76) 

where C is arbitrary. The form is the same as (J5U)) that we obtained in the model problem. 
We can also have reflections of the singularities at the intersection with the boundary. 
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The previous calculus are not valid if m — tanO and then A are equal to zero. In that case 
which is the static case, the propagation of singularities is rather different (see [8], [9]) and 
there is no reflection (cf. [7]). 

That study of the propagation of singularities has been done in a particular case of isotropic 
shell. The equations of vibration of a cylindrical shell have constant coefficients and then 
the propagations are along straight lines. The general case is more complicated: 
The segments which constitute the essential spectrum could be different in every point of the 
surface S. The characteristic curves are some pieces of curves and the propagations along 
them would depend on the reached point. For example, it is conceivable that the value of 
the spectral parameter which is given, will be go out of the essential spectrum in some point 
and that the propagation will stop. 

Moreover, if some propagation reaches a hyperbolic point at the edge of a shell with four 
directions of propagation, then we will have several possibilities for the reflection and we do 
not know what will happen. 
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